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BROAD LINE REGION OF FLAT SPECTRUM RADIO 
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ABSTRACT 

The broadband spectral energy distributions (SEDs) of 18 GeV flat spectrum 
radio quasars (FSRQs) are collected and compiled from literature, in which both 
the jet emission and the accretion disk radiation can be observed, in order to 
investigate the correlations among their jet power {Pjet), accretion disk luminosity 
(Ldisk), and luminosity of broad line region (BLR, Lblr)- On the basis of the SED 
hts with the jet radiation and accretion disk radiation models, we calculate Pjet 
and Ldisk- No correlation between Pjet with either Ldisk or Lblr is found. With 
a sub-sample of Lblr for 13 GeV FSRQs, it is observed that Lblr is strongly 
correlated with their Ldisk- We also study the BLR covering factors (LBLu/Ldisk) 
of the GeV FSRQs in our sample, averagely which are smaller than that of the 
large samples of radio-loud and radio-quiet quasars. Pjet of some GeV FSRQs 
is higher than Ldisk, but Pjet of all the GeV FSRQs is lower than the accretion 
power of black hole (BH), which is estimated by Mc^ = Ldisk/0.1, indicating that 
the total accretion power of BH is sufficient to drive the jets in these sources; 
however the uncorrelation between Ldisk and Pjet of the GeV FSRQs may suggest 
that their jets are launched by the Blandford—Znajek process via extracting the 
rotational energy of BH. Using the LBLR^Ldisk relation of the GeV FSRQs, we 
estimate Ldisk of a BL Lac sample with their Lblr- A comparison of Lblr and 
Eddington ratio (Ldisk/LEdd) among BL Lacs, very radio-loud NLSl galaxies. 
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and FSRQs is also presented. It is found that along with the BL Lac-NLSl- 
FSRQ sequence Fblr and L^isk/LEdd increase, which may correspond to the 
change of the accretion disk structure and the transformation of the dominant 
mechanism for jet launching. This is also consistent with the division of their 
parent populations, i.e., low-excitation radio galaxies and high-excitation radio 
galaxies. 

Subject headings: quasars: general—BL Lacertae objects: general—galaxies: 
Seyfert—galaxies: jets—accretion disks 


1. Introduction 

It is generally believed that there is a supermassive black hole (BH) at the center of an 
active galactic nucleus (AGN); the gravitational potential of the BH is the ultimate energy 
source of the AGN. The accreted matter of a BH forms an accretion disk and loses angular 
momentum through viscous or turbulent processes (Rees 1984). Broad emission lines are 
one of the dominant features of many AGN spectra, which are produced in a ‘broad line 
region’ (BLR) nearby the BH. The BLR is widely acknowledged to be photoionized by the 
continuum from the accretion disk (e.g., Kwan & Krolik 1981; Ferland et al. 1992). In order 
to explain the special radiation properties of BL Lac objects (BL Lacs), Blandford & Rees 
(1978) proposed that there are relativistic jet structures in these sources. The jet launching 
is reported to be connected with the central BH, accretion disk, and corona in a source 
(Armitage & Natarajan 1999; Merloni & Fabian 2002; Wang et al. 2004; Gao 2004, 2014; 
Yuan & Narayan 2014; Wilkins & Gallo 2015; Ghen et al. 2015). The jets may be powered via 
the Blandford—Payne (BP; Blandford & Payne 1982) process by releasing the gravitational 
energy of accreting matter that moves toward the BH, and/or via Blandford—Znajek (BZ; 
Blandford & Znajek 1977) mechanism by extracting the rotational energy of BH. 

So far, most conhrmed extragalactic GeV-TeV emitters are blazars, a sub-sample of 
radio-loud AGNs with relativistic jets. Blazars are divided into BL Lacs and flat spectrum 
radio quasars (FSRQs) according to their emission line features. Generally, the broadband 
spectral energy distributions (SEDs) of BL Lacs can be well explained with the synchrotron 
radiation and the synchrotron self-Gompton (SSG) scattering of the relativistic electrons in 
jets (Maraschi et al. 1992; Ghisellini et al. 1996, 2010; Zhang et al. 2012; Liao et al. 2014; 
Yan et al. 2014); however, FSRQs need a black body spectrum component and a component 
from inverse Gompton (IG) scattering of the external photons to ht their broadband SEDs 
(e.g., Sikora et al. 1994; Ghisellini et al. 1996, 2010; Zhang et al. 2014; Gao & Wang 2013; 
Ghen et al. 2010, 2012). Gompared with BL Lacs, FSRQs have prominent emission lines in 
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their spectra, and sometimes there are signihcant blue bumps at ~ 10^® Hz in their observed 
broadband SEDs, which is believed to be the thermal emission of accretion disk. However, 
often no signihcant blue bump is observed in the SEDs since the nonthermal radiation from 
jet overwhelms the thermal mission from accretion disk, especially when the sources are in 
the high state (e.g., for 3C 454.3, Jorstad et ah 2013). Hence, the GeV FSRQs in which the 
broadband SEDs are dominated by the jet emission and the signihcant blue bumps of the 
thermal emission from accretion disk are observed in their SEDs, are the good candidates 
for investigating the correlations among their jet power, accretion disk radiation, and BLR 
luminosity. 

The hux ratios of emission lines from BLR are suggested to be approximately constant 
among diherent AGNs (Francis et ah 1991; Vanden Berk et ah 2001). If setting the emission 
line luminosity of Lya to 100, the relative luminosities of H/9, Ha, Mg ii, and G iv lines are 
22, 77, 34, and 63, respectively, and then the total emission line luminosity is 555.8 (Francis 
et ah 1991; Gelotti et ah 1997). Therefore, the luminosity of BLR can be roughly estimated 
with the luminosities of one or some emission lines. Using the continuum luminosity and 
the emission line luminosity in Shen et al. (2011), Galderone et al. (2013) calculated Ldisk 
and Lblr, respectively, and reported that log(Ldisk/-hBLR) ~ 0.92 ± 0.28 — 1.08 ± 0.28, i.e., 
the covering factor of BLR is probably equal to 10 percent; therefore, the value of IOLblr 
is always used as a proxy to the accretion disk luminosity (Ldisk)- Sometimes O.lLdisk is also 
used to estimate the luminosity of BLR in blazars when authors have obtained Ldisk and 
want to calculate the energy density of BLR (e.g., Ghisellini et al. 2010; Zhang et al. 2014). 
However, so far it has not be examined whether the TBLu/Tdisk ~ 0.1 relation is still valid 
for blazars, especially for the GeV-TeV blazars. 

In this paper, a sample of SEDs for GeV FSRQs is compiled to study their jet and 
accretion disk relations and their covering factors of BLRs. The sample selection and data 
are presented in Section 2. The models and results of SED htting are described in Section 3. 
The study of the covering factors of BLRs is in Section 4. The possible connection between 
the jet and accretion disk is shown in Section 5. Further investigation to the correlations 
among the jet, BLR, and accretion disk is described in Section 6. A simple discussion on 
the anisotropic radiation of accretion disk, the jet formation mechanism for GeV AGNs, and 
the implications of unihcation for radio loud AGNs is given in Section 7. A summary for 
our results is in Section 8. Throughout, Hq=73 km s“^ Mpc“\ Da = 0.73 and flm = 0.27 
are adopted. 
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2. Sample and Data 
2.1. Sample Selection and SED Data 

We select the Fermi/LAT FSRQs with conhrmed redshift, simultaneously or quasi- 
simultaneously observed broadband SEDs, and obvious blue bumps in their SEDs. There 
are 18 sources (hereafter GeV-FSRQs) as shown in Figured! seven GeV-FSRQs (3G 273, PKS 
0420-01, PKS 0528+134, PKS 0727-11, IJy 1308+326, PKS 1510-089, PMN 2345-1555) 
from Zhang et ah (2014), eight GeV-FSRQs (PKS 0208-512, B3 0650+453, PKS 1127-145, 
PKS 1508-055, TXS 1846+322, TXS 2141+175, PKS 2144+092, PKS 2204-54) from Ghis- 
ellini et ah (2010), and three GeV-FSRQs from other literature (3G 454.3 in Vercellone et 
ah 2012; PKS 2123—463 in D’Ammando et ah 2012; PKS 2142—758 in Dutka et ah 2012). 


2.2. Luminosity and Energy Density of BLR 

The BLR luminosity (Lblr) is calculated with the method proposed by Gelotti et al. 
(1997) using the relative weights of different emission lines and the fluxes of emission lines 
whose observation data are available. The BLR luminosities of GeV-FSRQs in our sample are 
taken from Xiong & Zhang (2014), except for PKS 0528+134, for which Lblr is calculated 
using the flux of emission line G iv and Equation (1) in Gelotti et al. (1997). If there are 
more than one value available for Lblr in Xiong & Zhang (2014), the average of Lblr is 
used in this paper. A sub-sample of Lblr for 13 GeV-FSRQs in our sample is obtained, and 
the values of Lblr are shown in Table 1. We also consider the correction of Lblr since the 
different cosmological parameters are used in those works. 

The energy density of BLR at rest frame is calculated by Rbcr = 4 ;^^, where Rblr is 
the radius of BLR. The radius of BLR is estimated by the relations of log Rblr = “21.3l|g + 
0.519lg;Qgg logLsioo/lO^^^erg s“^ (Rblr in units of light days, Bentz et al. 2009) and Lh /3 = 
(1.425 ± 0.007) X 10^^( igUergVi srg s“^ (Greene & Ho 2005), where L5100 and 
Lh^ are the luminosities of continuum at 5100 A and H/5 line, respectively. If there is no 
observational flux of H/5 line available, Lh /3 is derived by the luminosities of other emission 
lines or Lblr using the relative flux ratio among different emission lines (Francis et al. 1991; 
Gelotti et al. 1997); please refer to Table 1 for more details. 
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2.3. Black Hole Mass 

The mass of BH, along with mass accretion rate, is a fundamental property of AGNs. 
Traditionally, the BH masses are estimated by assuming that the broad-line clouds are 
virialized (e.g., Krolik et ah 1991; Wandel et ah 1999; Wu et ah 2004; Greene & Ho 
2005). Another approach to estimate the mass of BH is to use the correlation between BH 
mass and stellar velocity dispersion by direct or indirect measurement of stellar velocity 
dispersion (Gebhardt et ah 2000; Ferrarese & Merritt 2000; Wu & Han 2001), which is 
always used to estimate the BH masses in BL Lacs (Woo & Urry 2002; Wu et ah 2002). 
The BH masses can also be derived using the relation between the BH mass and host-galaxy 
bulge luminosity (e.g., Laor 2001). There are 13 GeV-FSRQs in our sample for which the 
BH masses are obtained from literature, as listed in Table 2; the BH masses of 11 sources 
are estimated with the traditional virial method. For the other two soueces, one is derived 
with the relation of BH mass and bulge luminosity and the other one is estimated with the 
variability timescale. 


3. SED Modeling and Results 

As mentioned in Section 2, there are signihcant blue bump components at rvj 10^5 Hz 
in the observed SEDs of GeV-FSRQs in our sample, which are due to the thermal emission 
of accretion disk. Therefore, our model includes both the jet emission and accretion disk 
emission and the details of the model are described as follows. 


3.1. Accretion Disk Radiation 

The standard accretion disk spectrum in Davis & Laor (2011) is used to explain the 
signihcant bumps at ultraviolet band in the observed SEDs. The parameters include the 
inside and outside radius (Rin and Rout) of the accretion disk, BH mass, Eddington ratio, 
and inclination to the line of sight 9. The BH mass is hxed and the hxed values of BH masses 
for the sources are listed in Table 20. We take Rout = SOORg (Rg is the Schwarzschild radius) 
and cos 6 = 1, and then vary both Rin (from 0.5Rg to 4.5Rg) and the Eddington ratio to 
model the accretion disk emission at ultraviolet band. 


^The BH masses of the five sources for which no BH mass is found in literature are taken as the average 
of others, i.e., Mbh = 1.1 x 10® Mq. 



3.2. Jet Radiation 


The issue that the radiation regions are inside or outside the BLRs is very important to 
understand the properties of radiation regions, but highly debated. Considering the limits of 
Klein-Nashina (KN) effect to explain the observation data at TeV band and the absorption 
of 7 -rays via photon-photon pair production above 10 GeV (Liu & Bai 2006) of BLR photon 
held, some authors suggested that the emission regions should be outside the BLR and the 
photons from the torus dominate the IC process (Sikora et ah 2009; Tavecchio & Mazin 
2009; Tavecchio et ah 2013). However, Poutanen & Stern (2010) reported that the spectral 
breaks in the 2—10 GeV range can be well reproduced by the absorption of y-rays via photon- 
photon pair production on He ii Lyman recombination continuum and lines, implying that 
the radiation regions of blazars are inside the highest ionization region of BLR. Isler et al. 
(2013) also reported that the correlation of the increased emission line flux with millimeter 
core ejections and y-ray, optical, and ultraviolet flares implies that the BLR extends beyond 
the y-emitting region during the 2009 and 2010 flares for 3C 454.3. Sometimes the different 
timescales of variability are used to investigate the place of radiation region (e.g., Sikora et 
al. 2009). Some authors used two-zone model to explain the observed SED for the above 
controversies (e.g., Tavecchio et al. 2011). All the observed highest photon energies in our 
sample are smaller than 10 GeV, except for PKS 0727-11, which is ~20 GeV, and thus 
we still assume that the radiation regions are inside the BLRs and the external photon 
field is dominated by BLR, not torus. Therefore, the single-zone synchrotron-|-IC model is 
used to explain the jet emission, where the IC process includes both the SSC process and 
the external Compton scattering of the BLR photons (EC/BLR). Although the KN effect 
and the absorption of high energy gamma-ray photons by extragalactic background light 
(Franceschini et al. 2008) are not signihcant in the model calculation, we still take them into 
account. 

The radiation region is assumed to be a sphere with radius R, which is obtained with 
R = (5cAt/(l -|- z), where 6 is the beaming factor. At is the variability timescale and taken 
with At = 12 hr following our previous works (Zhang et al. 2012, 2014). The electron 
distribution is taken as a broken power law, which is characterized by an electron density 
parameter (A'"o, in units of cm“^), a break energy yb and indices (pi and P 2 ) in the range 
of ye ~ [ymin, ymax]- The indices of Pi and p 2 are derived from the spectral indices of the 
observed SEDs as reported by Zhang et al. (2012). y^ax can be only poorly constrained 
by the last data point in the SEDs and is usually dozens of times of yb, but it does not 
signihcantly affect our results. The energy density of BLR in the jet is boosted by a factor 
of P^ (P is the bulk Lorentz factor) and a corrected factor of 17/12 for Vblr proposed by 
Ghisellini & Madau (1996) is also considered because the BLR is not uniform and has an 
angle to the disk. So the derived energy density of the BLR photon fields in the comoving 
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frame is f/eLR = ^r^f/eLR- The values of ^7 /blr are listed in Table 1. For the sources 
whose T/blr are not available, the average of ^f/eLR = 2.73 x 10“^ erg cm“^ is taken. Since 
for blazars we are likely looking at the jet within the l/F cone, and that the probability is 
the highest at the rim of the cone, we take <5 = F in all the calculations, i.e., the viewing 
angle is equal to the opening angle of the jet. 

The free parameter set of our SED modeling is {S, <5, A^"o, 7b, 7min}, where B (in units 
of G) is the magnetic held strength of the radiation region. Following our previous works 
(Zhang et ah 2012, 2014), the minimization technique is also used to perform the SED 
hts. Since no errors are available for some data points in the selected SEDs, we estimate 
their errors with the average relative errors of the data points whose errors are available 
in order to obtain the goodness of the SED £t with the minimization technique. The 
average relative errors of the data for the bumps are derived separately (with a separation 
at z/ = 10^° Hz; e.g., Zhang et ah 2012, 2014) because the relative errors of the data for the 
synchrotron bump and SSC bump are usually much smaller than those for the EC bump. For 
3C 454.3, PKS 1510—089, and PKS 2123—463, no error of observation data below v = 10^° 
Hz is available, and we take 10% of the observation hux as the errors (Zhang et ah 2012; 
2014). For the details of this technique and htting strategies please refer to Zhang et al. 
(2012, 2014). 


3.3. Modeling Results 

Note that we only re-fit the jet emission for PKS 0528-1-134 and PMN 2345—1555 in 
this paper since the radiation of their accretion disks is taken as jet emission in Zhang et 
al. (2014), but the jet parameters for the other hve sources (3C 273, PKS 0420—01, PKS 
0727—11, IJy 1308-1-326, PKS 1510—089) in Zhang et al. (2014) do not change. The best 
SED hts are shown in Figure 1 and the derived model parameters in the la conhdence level 
are reported in Table 2. The distributions of jet htting parameters B, 6, Nq, 7 b, and 
are shown in Figure [2J 6 ranges from 7.4 to 20.6 with a median of ~14. B is from 1.3 to 11.5 
with a median of ~3.7. Most of 7min of the sources cluster within 10~100 with a median of 
~50; however, only one source (PKS 2142—758) needs a large 7min ~ 197 to ht the observed 
spectrum for the very hard spectrum at X-ray band. The distribution of 7b clusters within 
100~400, but the source TXS 2141-1-175 needs very small 7min and 7b to explain the observed 
spectrum at X-ray and GeV bands. The values of Nq mostly cluster in 3 x 10^ — 3 x 10'^ 
cm“^. From the htting spectra, we also calculate the luminosities of accretion disk radiation 
(Tdisk) and jet radiation (Tjet)- The values of Ldisk and Ljet are also reported in Table 2. 


4. Correlations Between Accretion Disk Luminosity and BLR Luminosity 


As described in Section 2.2, a sub-sample of BLR luminosity (Lblr) for the 13 GeV- 
FSRQs in our sample is obtained, and the luminosity of accretion disk (Ldisk) can be esti¬ 
mated from htting the blue bump. Lblr for the 13 GeV-FSRQs as a function of Ldisk is 
presented in Figure [31 Lblr is strongly correlated with Ldisk for the 13 GeV-FSRQs with a 
Pearson correlation coefficient of r = 0.85 and chance probability oi p = 2.7 x 10“^. The 
linear fit in the log scale gives logLBLR = (0.46 ± 6.04) -|- (0.96 ± 0.13) logLdisk- Here and 
hereafter, the bisector of the two ordinary least-squares lines of linear regression fit is used 
in order to avoid specifying independent and dependent variables (Isobe et ah 1990). With 
this correlation, one can estimate the accretion disk luminosity using the luminosity of BLR 
when the nonthermal emission of jet overwhelms the thermal emission of accretion disk. 

In order to compare the LBLR^Tdisk relation of GeV-FSRQs with the normal quasars, 
the sample of the 105,783 quasars from Shen et ah (2011) is compiled, which is taken in 
the Sloan Digital Sky Survey Data Release 7 quasar catalog. We only select the sources 
that have H/3 or/and Mg ii emission line^ with the signal to noise ratio of the lines larger 
than 10. To include radio properties, we compile two sub-samples with the radio loudness 
of Rrl = /e cm/f 2500 , where /e cm and /2500 are the flux density at rest-frame 6 cm and 2500 
A, and then we get 1999 very radio-loud quasars with Rrl > 100 (hereafter RLQs) and 937 
radio-quiet quasars with Rrl < 10 (hereafter RQQs). In Shen et ah (2011), the bolometric 
luminosities (Lboi) are derived from the continuum luminosities at 5100 A (L 5100 , 2 ; < 0.7), 
3000 A (L 3000 , 0.7 < 2 ; < 1.9), and 1350 A (L 1350 , 2 : > 1-9) with the correlations between Lboi 
and continuum luminosities in Richards et al. (2006). The values of Ldisk for the selected 
samples of RQQs and RLQs in Shen et al. (2011) are derived by ^Lboi (Galderone et al. 2013). 
The luminosities of BLRs for these quasars are also calculated with the fluxes of emission 
lines (H/3 or/and Mg ii) using the method proposed in Gelotti et al. (1997). The data of the 
selected RQQs and RLQs in Shen et al. (2011) are also presented in the TuLR^-hdisk plane, 
and Lblr are strongly correlated with Ldisk for the selected RQQs and RLQs, respectively. 
The linear fits in the log scale give logLBLR = —(2.19 ± 0.53) -|- (1.02 ± 0.01) logLdisk for 
RQQs and logLBLR = —(1.63 ± 0.59) -|- (1.01 ± 0.01) log Ldisk for RLQs, respectively. 

It can be found that the slope of the LBLR-Ldisk relation for GeV-FSRQs is consistent 
with that of both RQQs and RLQs, but the intercepts are slightly different among the 
three samples. We test whether the LsBR^Ldisk relations between GeV-FSRQs with RQQs 


^Compared with the emission line of C IV, the relative fluxes of H/3 and Mg II are more consistent with 
that reported by Francis et al. (1991), and the redshift range of the selected sample is also more consistent 
with ours. 



9 


and RLQs show any statistical difference with the two-dimensional Kolmogorov—Smirnov 
test (K-S test), which yields a chance probability of pxs- A two-dimensional K—S test 
probability larger than 0.2 (0.1 for one-dimensional) wonld strongly snggest no statistical 
difference between two samples, and a K—S test probability smaller than 10“^ would strongly 
suggest there is statistical difference between two samples. All the values of pxs derived with 
the two-dimensional (also the following one-dimensional) K—S tests between two samples are 
presented in Table 3. The results of K—S tests show that the 13 GeV-FSRQs in the Tblr^ 
-hdisk plane are not statistically distinct from the RQQs (pxs = 0.29), but they are marginally 
distinct from the RLQs (pxs = 0.045). 

To further investigate this issue, the distributions of Lblr and TdiskBas well as Tblr/ Ldisk 
for the 13 GeV-FSRQs and the selected samples of RQQs and RLQs are given in Figure ID 
Gonsidering the much different numbers of sources in the different samples, the distributions 
are normalized (i.e., the total number of the sources in the distribution has been set to unity) 
as shown in Figure 01 The mean of TBLu/Ldisk for the 13 GeV-FSRQs is 0.043, different from 
0.082 and 0.076 for RLQs and RQQs, respectively. The one-dimensional K-S tests indicate 
that the distribution of Lblr for the 13 GeV-FSRQs is analogous to that of both RQQs and 
RLQs, but the distribution of Ldisk for the 13 GeV-FSRQs is marginally different from that 
of both RQQs and RLQs in Shen et ah (2011), and thus the distribution of covering factor 
(LBLu/Ldisk) for the 13 GeV-FSRQs is different from that of the large quasar samples, either 
RQQs or RLQs. This also results in the similar slopes but the slightly different intercepts 
among the three samples in the TBLR^Ldisk plane. 


5. Correlations Between Accretion Disk Luminosity and Jet Power 

The jet power is estimated with the SED fitting parameters by assuming that they are 
carried by relativistic electrons, cold protons, magnetic fields, and radiation (e.g., Ghisellini 
et ah 2010; Zhang et al. 2012, 2013, 2014), i.e., Pjet = where P, = 2'kR^T‘^cU[ (the 

factor of 2 accounts for twin jets) and U[{i = e, p, P,r) are the powers and the energy 
densities associated with the emitting electrons, cold protons, magnetic fields, and radiation 
(e.g., Zhang et al. 2013, 2014), where the proton-electron pair assumption is used. The 
distributions of jet power and each ingredient power for the 18 GeV-FSRQs are shown in 
Figure O Pjet ranges within 5 x 10^^ — 8 x 10^® erg s“^. The distribution medians of jet 
ingredient powers are 5.85 x 10"^^ erg s“^, 9.14 x 10^® erg s“^, 2.31 x 10^^ erg s“^, and 


^The values of Ldisk only for the GeV-FSRQs in which Lblr is available are taken into account in this 
section, hence there are 13 sources, not 18. 
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3.08 X 10^® erg s“^ for P^, Pp, Pb, and Pr, respectively. The means of the jet radiation 
efficiency {Ej. = Pr/Pjet) and the jet magnetization parameter (a = Pb/(P p + Pe + Pr), Zhang 
et ah 2013, 2014) for the 18 GeV-FSRQs in our sample are 0.22 and 0.29, respectively, which 
is consistent with the result of our previous works that the FSRQ jets are highly magnetized 
with high radiation efficiency. 

The jet power (Pjet) as functions of the accretion disk luminosity (Tdisk) and the Ed¬ 
dington ratio (Tdisk/LEdd), and the jet radiation power (Pf) as a function of Tdisk are given 
in Figure O There are 10 sources in Zhang et al. (2014, hereafter FSRQs in Z14) that are 
not included in the selected sample due to lack of signihcant blue bump in their SEDs, for 
which the data of Lblr are available in literature, as listed in Table 4. There are also some 
GeV BL Lacs 0 in Zhang et al. (2012) for which Tblr is available in literature as listed in 
Table 4. Thus the data of the 10 FSRQs in Z14 and some BL Lacs in Zhang et al. (2012) 
are also shown in Figure |6l and Ldisk is derived from Lblr with the correlation between Ldisk 
and Lblr in Section 4 (Figure [3]). For comparison, the data of a GeV NLSl galaxy sample 
taken from Sun et al. (2015) are also given in Figure [6l Note that the values of Pjet and Pr 
for the 10 FSRQs in Z14, BL Lacs in Zhang et al. (2012), and GeV NLSl galaxies in Sun et 
al. (2015) have been multiplied by 2 to account for twin jets. 

In the Ldisk^Pjet plane, one can End that Ldisk is higher than Pjet for most of GeV 
FSRQs, but not so for BL Lacs in which the derived Ldisk is lower than Pjet. The data of 
GeV NLSl galaxies are in a broad range and are intermediate between GeV FSRQs and BL 
Lacs, as reported by Sun et al. (2015); hence the GeV NLSl galaxies are indeed a bridge 
between the FSRQs and BL Lacs. Assuming Ldisk = O.lMc^, where M is the accretion rate, 
the available accretion power of BH can be estimated by Mc^ = Ldisk/0.1. As shown in 
Figure |6](a), Pjet of all the FSRQs and GeV NLSl galaxies is lower than Mc^, but not so 
for most of BL Lacs. This may be due to that it is not a standard thin disk in the center 
of BL Lacs, but a radiatively inefficient disk (e.g., Zhang et al. 2014; Sun et al. 2015); in 
this case, Ldisk is lower than O.lMc^ (Yuan & Narayan 2014), i.e., Mc^ > Ldisk/0.1 for BL 
Lacs. Therefore, intrinsically Pjet of BL Lacs may be also lower than Mc^. No correlation 
between Ldisk and Pjet is observed for the 18 GeV-FSRQs in our sample, even taking the 10 
FSRQs in Z14 into account. However, if the data of BL Lacs and GeV NLSl galaxies are 
taken into account, Pjet would be strongly correlated with Ldisk with the Pearson correlation 
coefficient of r = 0.82 and chance probability of p = 1.5 x 10“^^. The linear fit in the log 
scale gives log Pjet = (18.82 ± 2.44) -|- (0.59 ± 0.05) logLdisk- Moreover, Pjet is also strongly 


"‘Including six sources, eight SEDs (two states with well constraints on ymin for BL Lacertae and W Com), 
they are Mkn 421, Mkn 501, lES 1218-1-304, W Com, PKS 2005—489, and BL Lacertae as given in Table 4 
and Zhang et al. (2012). 
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correlated with Ldisk/-hEdd for FSRQs, NLSl galaxies, and BL Lacs as shown in Fignre[6l^b); 
the Pearson correlation coefficient of r = 0.75 and chance probability of p = 4.7 x 10“®. 

In the Ldisk^-fr plane, Pj. for all the sonrces is lower than their Ldisk, except for one 
BL Lac. A tentative correlation between the two parameters seems to be observed for 
the 18 GeV-FSRQs in onr sample with the Pearson correlation coefficient of r = 0.48 and 
chance probability of p = 0.042. If taking the 10 FSRQs in Z14 into acconnt, the tentative 
correlation becomes slightly stronger with r = 0.50 and p = 0.007. It is interesting to find 
that BL Lacs, GeV NLSl galaxies, and FSRQs form a seqnence in the Ldisk^-G- plane; the 
linear £t in the log scale gives logPr = (6.16 ± 1.99) + (0.85 ± 0.04) logLdisk with a Pearson 
correlation coefficient of r = 0.93 and chance probability of p ~ 0. 

These resnlts may indicate that the total accretion power of BH is snfficient to lannch 
the jets in GeV-FSRQs, at least for the jet radiation in low state when the accretion disk 
radiation is not overwhelmed by the jet radiation; however the lack of strong correlation 
between Ldisk and Pjet niay imply that their energy sonrces are not exactly the same for GeV- 
FSRQs. Bnt the strong correlations between Pjet with Ldisk and Ldisk/-f"Edd for FSRQs, GeV 
NLSl galaxies, and BL Lacs may imply that the Eddington ratio may be the fnndamental 
in the nnihed framework among different types of GeV AGNs. 

Note that above resnlts are based on the traditional assnmption of one cold proton for 
one relativistic electron in the calcnlations of jet powers. There is no way to directly measnre 
the ratio of protons to electrons in the jets. The proton-electron pair assnmption is widely 
adopted in the calcnlations of jet powers in blazars (e.g., Ghisellini et ah 2010; Zhang et ah 
2012, 2014; Ghen et ah 2012). We check another extreme case that the jet power is carried 
by positron-electron pairs, magnetic field, and radiation, bnt no protons. In this scenario, 
the jet powers are still nncorrelated with Ldisk for the 28 FSRQs, and Pjet for the FSRQs 
and GeV NLSl galaxies are more lower than the accretion powers of BH, even for half of BL 
Lacs. 


6. Relations Among Jet Power, BLR Luminosity, and Accretion Disk 

Luminosity 

In order to fnrther stndy the relations among the jet, accretion disk, and BLR, Pjet 
and Pr as a fnnction of Lblr for the 13 GeV-FSRQs in onr sample are shown in Fignre 
[3 The data of the 10 FSRQs in Z14, some BL Lacs in Zhang et al. (2012), and the GeV 
NLSl galaxies in Snn et al. (2015) are also presented in Fignre [3 No correlation for FSRQs 
in the Pjet^-f^BLR plane is fonnd. Bnt there is a weak tentative correlation between P^ and 
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-^BLR for the 23 GeV FSRQs with the Pearson correlation coefficient of r = 0.48 and chance 
probability of p = 0.02. This tentative correlation may be due to that the emission of GeV 
FSRQs is dominated by EG process, for which the photon held is from BLR. R is similar to 
the Ldisk^-Pr plane, BL Lacs, GeV NLSl galaxies, and FSRQs also form a sequence in the 
L^ui-Pr plane. 

Then we make multiple linear regressions of the 13 GeV-FSRQs for Ldisk, -hsLR, -Pjet and 
for Ldisk, -hnLR, -Pr, respectively. The multiple linear regressions yield logLdisk = (—9.7 ± 
12.3) + (0.85 ± 0.16) logLBLR + (0.39 ± 0.24) logPjet with R-Square of 0.73 and logLdisk = 
(2.3 ±10.5) ± (0.83 ±0.20) logLBLR± (0.15 ±0.25) logPr with R-Square of 0.67, respectively. 
The predicted disk luminosity of Ldisk’ derived with the above htting functions of the multiple 
linear regressions, as a function of the observed accretion disk luminosity (Ldisk) is shown 
in Figure [Hi The best linear hts in the log scale give -f"disk^° with R-Square of 

0.75 and oc Ldisk^° with R-Square of 0.70 for using Pjet and Pr, respectively. These 
results further conhrm that Ldisk for the 13 GeV-FSRQs is tight correlated with Lblr, but 
with neither Pjet nor P^.. This indicates that Ldisk has the same energy origin as Lblr, but 
not Pjet; Ldisk is all from accretion process, but Pjet rnay be dominated by extracting the 
rotational energy of BH. Hence the dominant mechanism of jet launching in these GeV- 
FSRQs may be BZ process as discussed in our previous works (Zhang et al. 2014; Sun et ah 
2015). 


7. Discussion 

7.1. Anisotropic Radiation of Accretion Disk 

As described in Section 4, the different covering factor pBLu/Pdisk of BLR for the 13 
GeV-FSRQs with the selected RQQs and RLQs should be due to the different distributions of 
their Ldisk- In order to further compare the states of accretion disk among the three samples, 
the normalized distributions of BH masses (Mbh) and Eddington ratio (Ldisk/pEdd) are also 
shown in Figure |9l The results of K—S tests declare that the distribution of Mbh for the 
13 GeV-FSRQs is similar to that of RLQs and RQQs, but the distribution of Ldisk/pEdd for 
the 13 GeV-FSRQs is different from both RLQs and RQQs. Therefore, the main reason for 
the different TBLn/Pdisk and Ldisk/pEdd of the 13 GeV-FSRQs with RLQs and RQQs is the 
different distributions of their Ldisk- 

Note that the GeV-FSRQs in our sample belong to blazars, in which the relativistic 
jets (also the normal of the accretion disk) are thought to be close to the line of sight. 
With the same intrinsic luminosity, one may observe the higher flux of accretion disk for the 
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GeV-FSRQs than other sources if the radiation of the accretion disk is anisotropic (Liu Sz 
Zhang 2011). Assuming that the emitting flux of accretion disk depends on the inclination 
angle, i.e., F = FqCos^, where 6 is the angle between the line of sight and the normal of 
the accretion disk, the inferred intrinsic luminosity of accretion disk should be 
(see details in the Appendix), i.e., = Ldisk/2 for the GeV-FSRQs. Assuming that the 

inclination of accretion disk for RLQs and RQQs is the average value of ~ 60°, the derived 
isotropic luminosity of accretion disk is just the inferred intrinsic luminosity of accretion disk 
for RLQs and RQQs, i.e., = Ldisk- One would observe the similar distributions of BLR 

covering factor {Lblr/L^ isk) ^^e GeV-FSRQs with RLQs and RQQs if considering the 
inclination effect and anisotropic radiation of accretion disk; the BLR covering factor of the 
GeV-FSRQs would be close to 0.1 with a mean of 0.086, similar to RLQs and RQQs in Shen 
et ah (2011). 

Recently, Richards & Lister (2015) reported that the jets of radio-loud (RL) NLSl 
galaxies are aligned at moderately small angles to the line of sight, which is similar to 
blazars. Hence the intrinsic luminosity of accretion disk for GeV NLSl galaxies should 
have the same correction factor as FSRQs and BL Lacs, i.e., = Ldisk/2. Taking the 

anisotropic radiation of accretion disk into account in the L^^^y-Pjet and planes, the 

results are given in Figure ITOl For most of the sources, Pjet is higher than Ldiski but most of 
the FSRQs and NLSl galaxies still have higher accretion power of BH (Me ^ = ^Lk/O-l) than 
Pjet, indicating that the pure accretion power is still enough to drive their jet. In the L^^^y-P,. 
plane, the sources still form a sequence and are also closer to the equality line, implying that 
the radiation power of jet is already comparable with the luminosity of accretion disk for the 
GeV sources. 

Following our previous works, the SED hts in this paper are still under the assumption 
of (5 = F, i.e., the viewing angle is equal to the opening angle (1/F) of a jet, since the 
probability is the highest when looking at the jet at the angle of 1/F as described in Section 
3.2. It is well known that there is Doppler boosting effect in the radiation of balzars. If the 
viewing angle is larger than the opening angle of the jet, the leptonic models would not be 
able to explain the observation data. In another scenario, 5 = A * F when the viewing angle 
is smaller than the opening angle of the jet, where the coefficient of A ranges from 1 to 2, the 
calculated jet powers with the derived parameters would be lower than the previous results 
(under the assumption of 6 = F). That means the jet powers may be overestimated under 
the assumption of <5 = F to fit the SEDs; under this case, the jet powers are further lower 
than those that derived with the assumption of 5 = F, strengthening the above conclusion 
that the jet powers are lower than the accretion power of BH. 
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7.2. Eddington Ratio and Jet Formation of GeV AGNs 

The radiation of FSRQs and BL Lacs as well as the conhrmed GeV NLSl galaxies (Abdo 
et ah 2009; Paliya et al. 2013; 2014; Sun et ah 2015; Angelakis et al. 2015) is believed to 
be dominated by the emission of relativistic jets. The observed spectral properties of the 
GeV NLSl galaxies are analogous to GeV FSRQs with signihcant broad emission lines and 
the blue bump from accretion disk radiation. BL Lacs are dehned as the equivalent width 
(EW) of emission lines in the rest frame being less than 5 A, and no signihcant blue bump 
at ~ 10^^ Hz is observed in their SEDs. In recently years, some sources that are classihed as 
BL Lacs show some properties similar to FSRQs (Sbarufatti et al. 2005; Raiteri et al. 2007; 
Ghisellini et al. 2011; Giommi et al. 2012), implying that their thermal emission from the 
accretion disk may be overwhelmed by the jet emisison. Hence the emission line luminosity 
of the BL Lacs can be used as a proxy to investigate the accretion disk luminosity and then 
the Eddington ratio of these sources. We collect and compile a sample of Lblr for BL Lacs 
from literature, which includes 16 BL Lacs with available BH masses as listed in Table 4. 
Their Ldisk are also estimated by Lblr with the correlation between Ldisk and Lblr of the 
13 GeV-FSRQs in Figure [3l 

For NLSl galaxies, besides GeV NLSl galaxies PMN J0948+0022, PKS 1502+036, and 
IH 0323+342, in which their Ldisk are available and derived by htting the blue bump in 
SEDs with the standard accretion disk spectrum (Sun et al. 2015), the other very RL NLSl 
galaxies in Yuan et al. (2008) are also taken into account. Using the fluxes of emission lines 
reported in Yuan et al. (2008), Lblr of these NLSl galaxie^ are derived with the method 
proposed in Gelotti et al. (1997). The values of Ldisk/-f^BLR for the GeV NLSl galaxies 
PMN J0948+0022, IH 0323+342, and PKS 1502+036 are ~86, ~14, ~28, respectively. So 
40Lblr is used to estimate Ldisk of other RL NLSl galaxies in Yuan et al. (2008). It was 
suggested that the BH masses of NLSl galaxies derived from the data of the H/3 line may 
be underestimated (Gollin et ah 2006; Zhu et al. 2009; Galderone et al. 2013). Thus the 
BH masses of RL NLSl galaxies in Galderone et al. (2013) are used to calculate Eddington 
ratio. There are 17 very RL NLSl galaxies with available BH masses in Galderone et al. 
(2013). Adding IH 0323+342 (Sun et al. 2015), there are 18 RL NLSl galaxies as listed in 
Table 4. 

The distributions of Lblr and LdiskZ-f^Edd for 16 BL Lacs, 18 RL NLSl galaxies, and 
23 GeV FSRQs (including the 10 FSRQs in Z14) are shown in Figure [HI It is found that 
the distributions of Lblr form a sequence along BL Lac-NLSl-FSRQ, and Lblr of BL Lacs 
are intrinsically weak compared with that of RL NLSl galaxies and GeV FSRQs. It may 


®For IH 0323+342, the flux of emission line H/3 is taken from Yao et al. (2015). 
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be due to the intrinsically weak radiation of the accretion disk in BL Lacs, and thus it is 
easily overwhelmed by the jet emission. There is also a BL Lac-NLSl-FSRQ sequence in 
the Ldisk/-hEdd distribution. It can be found that Ldisk/-hEdd for most of the BL Lacs are 
lower than 0.01, but Ldisk/-hEdd for the GeV FSRQs is higher than 0.01 (see also Ghisellini 
et al. 2010; a division to separate BL Lacs and FSRQs). The Ldisk/FEdd distribution of 
NLSl galaxies is intermediate between FSRQs and BL Lacs, and is more analogous to that 
of FSRQs. 

As shown in Figured Tjet is strongly correlated with Ldisk and Ldisk/-^Edd if taking all the 
data of FSRQs, GeV NLSl galaxies, and BL Lacs into account, indicating the accretion rate 
is very important to drive jets. Zhang (2013) and Sun et al. (2015) also proposed that the 
accretion rate (see also Shen & Ho 2014) may be the fundamental in the unihed framework 
among different types of GeV AGNs; along with the BL Lac-NLSl-FSRQ sequence, the 
increase of the Eddington ratio corresponds to the change of the accretion disk structure, 
and then may correspond to the transformation of the dominant mechanism for jet launching. 
The jets of GeV FSRQs may be powered via extracting the rotational energy of BH with BZ 
process, which is also consistent with the uncorrelation between Ldisk and Pjet of the GeV 
FSRQs as described in Section 5. However, the jets of BL Lacs and NLSl galaxies may be 
produced via the BP and/or BZ mechanisms, depending on the structure and accretion rate 
of their accretion disks. Note that the above discussion does not consider the correction 
of anisotropic radiation of accretion disk mentioned in Section 7.1, but this effect does not 
change the results since the jets of GeV NLSl galaxies have the similar properties to that of 
blazars as discussed in Section 7.1. 


7.3. Implications of unification for RL AGNs 

The topic of the differences and unihcation among the different types of AGNs has been 
extensively covered and discussed. According to the unihed models for RL AGNs, BL Lacs 
are associated with center-brightened (FR I) radio galaxies, whereas FSRQs are usually fond 
in edge-brightened (FR H) radio galaxies (Urry & Padovani 1995). They may also correspond 
to the low-excitation radio galaxies (LERGs) and high-excitation radio galaxies (HERGs), 
although the low/high-excitation radio galaxies do not have one-to-one correspondence with 
the FR I—FR H categories (Hine & Longair 1979; Laing et al. 1994; Hardcastle et al. 2009). 
This is also consistent with numerous BL Lacs with FR H radio structure (e.g. Rector & 
Stocke 2001) and some quasars with FR I characteristic (e.g., Blundell & Rawlings 2001). 
Therefore, intrinsically LERGs and HERGs may be the parent populations of BL Lacs and 
FSRQs. Recently, some authors suggested that the intrinsic differences between HERGs 


and LERGs are due to the different modes of accretion, i.e., owing to the different accretion 
rates (Mingo et ah 2014; Fernandes et ah 2015); HERGs correspond to the radiatively 
efficient ‘cold mode’ accretion with higher accretion rates, whereas LERGs are the radiatively 
inefficient ‘hot mode’ accretion with lower accretion rates. The switch lies approximately 
at the Eddington rates 0.01~0.04, but the jet production is not switched off. However, the 
division between the two populations is not very clear, there is still some overlap with wider 
range of accretion rates for LERGs (Mingo et al. 2014). These results agree well with ours. 

As described above, some of the radio properties between HERGs with LERGs do not 
show clear differences and there is no dependence between jet power with radiative luminosity 
for HERGs, presumably implying that the process of power conversion into jet power is not 
solely controlled by the accretion rate; an additional process must be influencing the jet 
formation, such as the spin of BH (Mingo et al. 2014; Fernandes et al. 2015). This is also 
consistent with our results that there is no correlation between jet power with Eddington 
ratio for single population of FSRQs, and the extraction of rotation energy of BH may be 
important to power their jets. 

The large scale kinetic powers of RGs usually are estimated with the radio luminosity at 
151 MHz (e.g., Godfrey & Shabala 2013; Mingo et al. 2014) and they reflect the long-term 
average power injection from their jets over the whole RL lifetimes of the sources. Averagely, 
the large kinetic powers of LERGs are lower than that of HERGs (e.g., Mingo et al. 2014). 
Interestingly, the large scale kinetic powers of blazars is strongly correlated with their jet 
powers (Zhang et al. 2014), which reflect essentially their instantaneous activities at the 
present epoch, and the large scale kinetic powers of BL Lacs are lower than that of FSRQs 
(Zhang et al. 2014; Sun et al. 2015). The statistical associations of BL Lacs with LERGs 
and FSRQs with HERGs suggest that the current core activities of LERGs also have on 
average lower power outputs than HERGs. 

There are two possible scenarios for the above main differences between the two classes 
of sources, i.e., BL Lacs/LERGs and FSRQs/HERGs. One scenario is that the two classes 
are essentially the same sources in nature, but their accretion disks may stay in different 
accretion states at different times and thus produce different jet/core power outputs. In other 
words, the two classes of objects would be converted from each other when the accretion 
states change. This scenario agrees with the fact that their host galaxies seem to show 
similar properties (Urry & Padovani 1995; Seymour et al. 2007; Fernandes et al. 2015). 
However, a main consequence of this scenario is that their long-term average jet/core power 
outputs should be essentially the same, which contradicts with the observed higher large- 
scale jet kinetic powers of FSRQs/HERGs than BL Lacs/LERGs. The alternative is that 
the two classes are essentially different objects with BL Lacs/LERGs tending to produce 
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lower jet/core powers (i.e. with lower accretion rate and/or black hole spin), resulting in 
the lower large-scale jet kinetic powers in BL Lacs/LERGs. The apparent contradiction in 
this scenario with the indistinguishable host galaxies of these two classes may be resolved 
if they have different circum-nuclear environments resulting in different accretion supply on 
the average. Then the accretion disk instabilities and/or the different jet orientations may 
result in large scatter of the observed jet/core outputs at different times, in agreement with 
observations (e.g., Mingo et al. 2014). Future higher quality observations may test if their 
circum-nuclear environments are indeed different intrinsically. 


8. Summary 

A SED sample of 18 GeV-FSRQs, in which the signihcant blue bumps from the accretion 
disk radiation are observed and the BLR luminosities (Tblr) of 13 GeV-FSRQs are available, 
is collected and compiled from literature. On the basis of the SED hts with the jet radiation 
and accretion disk radiation models, we calculate the jet power (Tjet) and the accretion disk 
luminosity (Tdisk), investigate the correlations of the jet power and disk luminosity as well 
as the BLR luminosity, and also make a comparison of the BLR covering factor (TeLR/Tdisk) 
with the large RLQ and RQQ samples from Shen et al. (2011). Our results are summarized 
below. 

• -hRLR of the 13 GeV-FSRQs is strongly correlated with Ldisk, and thus Lblr can be 
used to estimate Ldisk when the radiation of accretion disk is overwhelmed by the jet 
emission. The mean of TBLu/Tdisk for the 13 GeV-FSRQs is 0.043, which is smaller 
than that of the RLQs (~0.082) and RQQs (~0.076) in Shen et al. (2011). Gonsid- 
ering the inclination effect and anisotropic radiation of accretion disk, the TBLu/Tdisk 
distribution of the 13 GeV-FSRQs would be same as that of RLQs and RQQs in Shen 
et al. (2011). 

• No correlation between Pjet with Ldisk and Lblr is found for the GeV FSRQs. Pjet 
for most of the GeV FSRQs is lower than Ldisk and all of them are lower than the 
accretion power of BH (Mc^ = Ldisk/0.1), indicating that the total accretion power of 
BH is sufficient to drive the jets in these sources, but the uncorrelation between Ldisk 
and Pjet of the GeV FSRQs in our sample may suggest that their jets are launched by 
the BZ process via extracting the rotational energy of BH. 

• Using the LBLR^Tdisk relation of the 13 GeV-FSRQs, we estimate Ldisk of a BL Lac 
sample with their Lblr; a very RL NLSl galaxy sample from Yuan et al. (2008) is also 
taken into account. It is found that along with the BL Lac-NLSl-FSRQ sequence Lblr 
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and the Eddington ratio (Edisk/EEdd) increase, which may correspond to the change 
of the accretion disk structnre and the transformation of the dominant mechanism for 
jet lannching. This is also consistent with the division of their parent populations, i.e., 
BL Lacs/LERGs and FSRQs/HERGs. 
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9. Appendix 

As shown in Figure [T21 the observed flux of accretion disk depends on the inclination 
angle, i.e., F = Fq cosd, where 9 is the angle between the line of sight and the normal of the 
accretion disk. The intrinsic luminosity of accretion disk is 

-^disk = Fq cos 9 sin 9d9 = 271 DIFq, (1) 

Jo Jo 

where Hl is the luminosity distance of the source. The derived isotropic luminosity is 

^disk = 4:71 DlF = iTiDlFo cos 9. (2) 

Therefore, the inferred intrinsic luminosity should have a correction factor with the derived 
isotropic luminosity, i.e., 

jB _ Fdisk 
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Table 1. The data for BLRs of 18 GeV-FSRQs in our sample 


Source 

FIuxh/3’^ 

Line’’ 

log Tblr’^ 

7?blr‘^ 

Fblr'^ 

3C 273 

1548 

H/5C97 

45.50 

238.3 

3.13 

3C 454.3 

43 

H/5C97 

45.57 

287.0 

2.51 

PKS 0208-512 

12 

Mg 

45.17 

190.9 

2.27 

PKS 0420-01 

8 

Mg 

44.88 

142.9 

2.07 

PKS 0528+134 

1.5 

C lyP'' 

45.06 

168.3 

2.28 

B3 0650+453 



44.24 

70.6 

1.95 

PKS 0727-11 






PKS 1127-145 

19.36 

Lya®"'® 

45.75 

287.4 

3.82 

IJy 1308+326 

3.98 

H/5C09 

44.99 

114.4 

4.20 

PKS 1508-055 



45.50 

267.3 

2.48 

PKS 1510-089 

53 

H/3C97 

44.72 

119.9 

2.04 

TXS 1846+322 



44.56 

98.7 

2.07 

PKS 2123-463 






TXS 2141+175 

20.95 

C IvO®4 

44.23 

45.0 

4.66 

PKS 2142-758 






PKS 2144+092 






PKS 2204-54 






PMN 2345-1555 



44.33 

78.1 

1.98 


'^The flux of emission line H/3 is in units of 10“^^ erg s“^ cm“^. Lblr 
is in units of erg s“^. The luminosities of line H/3 for B3 0650+453, PKS 
1508—055, and PMN 2345—1555 are derived with Lblr when it is used to 
estimate i?BLR- 

'’The line used to estimate the flux of emission line H/3 with the relative 
flux ratios. The superscripts denote the references-C97: Celotti et al. (1997); 
Pll: Palma et al. (2011); T12: Tang et al. (2012); C09: Chen et al. (2009); 
094: Osmer et al. (1994). 

‘^7?blr and I/blr are in units of light days and 10“^ erg cm“^, respectively. 







Table 2. Derived Parameters from Our SED Fits with the Single-zone Peptonic Model 


Source 

z 

Pi 

P 2 

5 

B 

[G] 

Tmin 

7b 

7max 

No 

[xlO^ cm~^] 

log ^jet^ 
[erg s“^] 

log idisk 

[erg s-1] 

log -Pjet 

[erg s-1] 

Mbh*^ 

[log Mq] 

3C 273 

0.158 

1.4 

4.0 

7.4+0.9 

8.5+1.6 


328+79 

2000 

5.6+5.3 

47.11+0.03 

46.92+0.10 

45.88+0.18 


3C 454.3 

0.859 

1.2 

3.46 

15.6+0.6 

5.1+0.8 

40+6 

137+32 

1800 

3.4+2.0 

48.49+0.02 

46.97+0.09 

46.50+0.13 

9 I7W02 

PKS 0208-512 

1.003 

1.3 

3.56 

15.8+0.7 

3.4+1.2 

40+13 

105+40 

3300 

6.9+6.9 

48.24+0.02 

46.46 

46.39+0.32 

g 21^04 

PKS 0420-01 

0.916 

2.64 

4.06 

12.8+0.7 

8.2+1.1 

51+11 

385+118 

5000 

2512+1388 

48.01+0.02 

46.18 

46.28+0.13 

g O3WO2 

PKS 0528+134 

2.07 

1.2 

4.2 

17.4+0.9 

2.9+1.1 

111+6 

230+86 

7000 

16.6+16.6 

49.01+0.03 

46.93 

46.66+0.22 

9 4X04 

B3 0650+453 

0.933 

1.28 

3.6 

14.1 + 1.0 

1.3+0.3 

95+16 

111+37 

4000 

15.8+15.8 

48.08+0.02 

46.02+0.17 

46.16+0.25 

31^312 

PKS 0727-11 

1.589 

1.8 

3.9 

20.6+1.2 

5.4+1.1 

50+12 

254+61 

5000 

15.8+10.6 

48.91+0.03 

46.59 

46.46+0.11 


PKS 1127-145 

1.184 

1.2 

4.56 

13.1+0.8 

11.5+1.2 

49+4 

213+32 

4000 

3.7+2.1 

48.15+0.03 

47.15+0.08 

46.35+0.09 

9.18^1® 

IJy 1308+326 

0.997 

2.2 

3.4 

12.6+0.9 

3.4+0.9 

44+10 

353+129 

8000 

490+451 

48.20+0.02 

46.09 

46.29+0.24 

g94009 

PKS 1508-055 

1.185 

1.2 

3.8 

17.0+1.1 

7.0+1.5 

32+14 

141+41 

5000 

0.9+0.7 

48.00+0.03 

47.04+0.17 

46.27+0.16 

gg^Z 09 

PKS 1510-089 

0.36 

1.12 

3.72 

11.0+0.5 

3.1+0.5 

20+12 

305+32 

1500 

0.4+0.2 

47.46+0.03 

45.77 

45.72+0.09 

g 65WO2 

TXS 1846+322 

0.798 

1.4 

3.76 

13.1+0.6 

2.5+0.7 

69+8 

206+72 

5000 

4.9+4.4 

47.63+0.03 

46.42+0.18 

45.82+0.22 

8 21S12 

PKS 2123-463 

1.67 

1.3 

3.8 

17.9+0.6 

3.6+0.6 

88+5 

243+56 

5000 

4.8+2.8 

48.82+0.02 

46.65 

46.37+0.09 


TXS 2141+175 

0.213 

1.4 

4.02 

10.3+0.6 

5.1+0.8 

2+1 

42+6 

8000 

15+8.1 

47.01+0.02 

46.05+0.03 

46.87+0.17 

8.98^°® 

PKS 2142-758 

1.139 

1.6 

3.6 

9.5+1.4 

1.7+1.0 

197+17 

422+188 

3000 

856+856 

48.56+0.03 

46.68 

46.67+0.26 


PKS 2144+092 

1.113 

1.4 

4.0 

14.3+1.0 

3.8+1.2 

74+8 

175+66 

4000 

11.5+11.5 

48.04+0.03 

46.29 

46.12+0.24 


PKS 2204-54 

1.215 

1.4 

4.4 

14.4+0.9 

5.7+1.2 

52+5 

205+46 

4000 

6.5+4.4 

47.95+0.03 

46.52 

46.12+0.14 


PMN 2345-1555 

0.621 

2.0 

4.0 

13.8+1.0 

2.6+0.8 

46+11 

141+55 

8000 

9.3+9.3 

47.50+0.04 

45.42 

45.83+0.25 

8.16®!® 


‘^The errors of jet radiation luminosity are estimated with the errors of the two peak luminosities only. 

masses of the sources. The superscripts denote the references—F04: Fan Sc Cao (2004); W02: Woo Sc Urry (2002); X04: Xie et al. 2004; S12: Shaw et al. (2012); Z09: Zhou 
Sc Cao; C09: Chen et al. (2009); and F03: Falomo et al. (2003). 
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Table 3. K—S Test Probability of for fhe Parameter Distributions among the 13 

GeV-FSRQs, RLQs, and RQQs 


Parameters'^ 

FSRQs & RLQs 

FSRQs & RQQs 

RLQs & RQQs 

l^disk'l^BLR 

0.045 

0.29 

1.3E-24 

Lbim 

0.38 

0.59 

4.9E-12 

-^disk 

0.06 

0.07 

7.1E-7 

l^BLR/idisk 

2.3E-4 

0.002 

3.9E-11 

Mbh 

0.13 

0.45 

1.5E-27 

idisk/l^Edd 

2.0E-5 

0.015 

6E-37 


“pKS larger than 0.1 (0.2 for two-dimensional K—S test) would strongly suggest no statistical difference 
between two samples, whereas pKS smaller than 10“^ would suggest there is significantly statistical difference 
between two samples. 

^The two-dimensional K—S test is only used for the Ldisk“-bBLR relation, and the one-dimensional K—S tests 
are given for the other parameter distributions between two samples. 







Table 4. BLR Luminosities and BH Masses of RL AGNs 



RL NLSl galaxies^ 



BL Lacs^ 


GeV FSRQf 

s in Z14‘= 


Source 

z 

log LblR 
[erg s”^] 

ATbh 

[log Mq] 

Source 

z 

log LblR 
[erg s“^] 

ATrh 

[log Mq] 

Source 

z 

log LblR 
[erg s”^] 

Mbh 

[log Mq] 

1H0323+342 

0.0629 

43.53 

8.6 

PKS 0521-36 

0.055 

42.68 

8.6 

S4 0133+47 

0.859 

44.42 

8.3 

J0814+5609 

0.509 

44.28 

8.4 

PKS 0829+46 

0.174 

42.57 

8.68 

4C 28.07 

1.213 

45.30 

9.22 

J0850+4626 

0.523 

43.93 

8.8 

PKS 0851+202 

0.306 

42.83 

8.86 

PKS 0454-234 

1.003 

44.39 

9.17 

J0902+0443 

0.532 

44.09 

8.5 

TXS 0954+658 

0.367 

42.45 

8.53 

S4 0917+44 

2.19 

45.77 

9.29 

J0948+0022 

0.584 

44.03 

9.1 

PMN 1012+063 

0.727 

42.89 

8.5 

4C 29.45 

0.729 

44.66 

8.61 

J0953+2836 

0.657 

44.01 

8.3 

PKS 1057-79 

0.581 

43.76 

8.8 

3C 279 

0.536 

44.48 

8.28 

J1037+0036 

0.595 

43.78 

9.1 

Mkn 421 

0.031 

41.69 

8.67 

PKS 1502+106 

1.839 

45.23 

8.98 

J1146+3236 

0.465 

44.11 

8.8 

lES 1218+304 

0.182 

42.06 

8.58 

B2 1520+31 

1.487 

44.88 

8.92 

J1238+3942 

0.622 

43.76 

8.2 

W Com 

0.102 

42.14 


4C 66.20 

0.657 

44.40 

9.14 

J1246+0238 

0.362 

43.77 

8.7 

PKS 1519-273 

1.294 

43.53 

8.8 

PKS 2325+093 

1.843 

45.19 

8.7 

J1305+5116 

0.785 

45.17 

9.2 

Mkn 501 

0.034 

42.20 

9.03 





J1435+3131 

0.501 

43.98 

8.8 

PKS 1749+096 

0.322 

43.70 

8.66 





J1443+4725 

0.703 

44.30 

8.4 

S5 1803+78 

0.68 

44.85 

8.36 





J1505+0326 

0.408 

43.39 

8.3 

TXS 1807+698 

0.051 

42.00 

8.7 





J1548+3511 

0.478 

44.26 

8.7 

PKS 2005-489 

0.071 

42.04 

8.76 





J1634+4809 

0.494 

43.86 

8.9 

BL Lacertae 

0.069 

42.52 

8.61 





J1644+2619 

0.144 

43.14 

8.2 

PKS 2240-260 

0.774 

43.46 

8.6 





J1722+5654 

0.425 

44.06 

8.1 










to 


is estimated with the fluxes of emission lines in Yuan et al. (2008). The BH masses of these sources are taken from Calderone et al. (2013). 
The data are taken from Ghisellini et al. (2011) and Sbarrato et al. (2012). 

The values of and are taken from Xiong & Zhang (2014). have been corrected for the different cosmological parameters. 
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Fig. 1.— Observed SEDs (black data points) with the model best hts (lines) for the 18 
GeV-FSRQs in our sample. The black solid lines are the sum of emission from the jet and 
the accretion disk (red dashed lines). 



Fig. 2.— Distributions of the magnetic held strength {B), the beaming factor (5), the 
minimum Lorentz factor of electrons (7min), the break Lorentz factor of electrons (7b), and 
the electron density paramete {Nq) for the 18 GeV-FSRQs in our sample. 
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Ldisk [erg 

Fig. 3.— Fblr as a function of L^isk for the 13 GeV-FSRQs (black solid circles) in our 
sample. The data of the RLQs (blue opened triangles) and RQQs (red opened squares) are 
from Shen et ah (2011). The lines are the bisectors of the two ordinary least-squares lines of 
linear regression hts to the data for GeV-FSRQs (black solid line), RLQs (blue dashed line), 
and RQQs (red dash-dotted line), respectively. 



Lblr [erg s'] [erg s’’] 

Fig. 4.— Normalized distributions of the BLR luminosity (Lblr), the disk luminosity (Ldisk), 
and the BLR covering factor {L^lr/L^ isk) for the 13 GeV-FSRQs (black solid lines) in our 
sample. The data of the RLQs (blue dashed lines) and RQQs (red dotted lines) are from 
Shen et ah (2011). 

















































30 



Pje, [erg s '] [erg s '] [erg s '] [erg s '] P, [erg s '] 

Fig. 5.— Distributions of the jet power (Pjet) and the powers associated with relativistic 
electrons {Pe), cold protons (Dp), Poynting flux (Db), radiation component (Pr) for the 18 
GeV-FSRQs in our sample. 





Fig. 6.— Pjet as functions of Ldisk (Panel (a)) and Ldisk/pEdd (Panel (b)), and Pr as a function 
of Fdisk (Panel (c)) for the GeV-FSRQs (black solid circles) in our sample. The data of the 
10 GeV FSRQs (black opened circles) in Zhang et ah (2014), the BL Lacs (blue triangles) 
with available Lblr in Zhang et ah (2012), and the GeV NLSl galaxies (magenta stars) 
from Sun et ah (2015) are also presented for comparison. The solid lines are the bisectors of 
the two ordinary least-squares lines of linear regression hts for all the sources in the hgures, 
including FSRQs, NLSl galaxies, and BL Lacs. 
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Lblr [erg s'] 


Fig. 7.— Pjet and Pj. as a function of Fblr for the 13 GeV-FSRQs (solid circles) in our 
sample. The data of the 10 GeV FSRQs (black opened circles) in Zhang et al. (2014), the 
BL Lacs (blue triangles) with available Lblr in Zhang et al. (2012), and the GeV NLSl 
galaxies (magenta stars) from Sun et al. (2015) are also presented for comparison. 



Fig. 8.— The relations of the observed accretion disk luminosity (Ldisk) and the predicted 
luminosity (T)((g]^) that are derived with the functions of the three element linear regression of 
logLj^C®^ = -9.7 +0.85 log Lblr +0.39 log Pjet and = 2.3 + 0.83 log Lblr + 0.15 logPr, 

respectively. The black solid line and the red dashed line are the best hts for Ldisk with 
and respectively. 
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Fig. 9.— Normalized distributions of BH masses and Ldisk/FEdd for the 13 GeV-FSRQs 
(black solid lines) in our sample. The data of the RLQs (blue dashed lines) and RQQs (red 
dotted lines) are from Shen et ah (2011). 




Fig. 10.— Pjet and Pj. as a function of the anisotropic luminosity of accretion disk (L^j^j^) 
for the 18 GeV-FSRQs (black solid circles) in our sample. The data of the 10 GeV FSRQs 
(black opened circles) in Zhang et ah (2014), the BL Lacs (blue triangles) with available 
-f^BLR in Zhang et al. (2012), and the GeV NLSl galaxies (magenta stars) from Sun et al. 
(2015) are also presented for comparison. 
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Fig. 11.— Distributions of L^lr and L^isk/LEdd for BL Lacs, very RL NLSl galaxies in Yuan 
et al. (2008), the GeV FSRQs in our sample and in Zhang et ah (2014). 
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Fig. 12.— Illustration of anisotropic radiation for accretion disk. 






































